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Abstract 

Several aspects of mathematical astrobiology are discussed. It is ar- 
gued that around the time of the origin of life the handedness of biomolecules 
must have established itself through an instability. Possible pathways of 
producing a certain handedness include mechanisms involving either au- 
tocatalysis or, alternatively, epimerization as governing effects. Concepts 
for establishing hereditary information are discussed in terms of the the- 
ory of hypercycles. Instabilities toward parasites and possible remedies 
by invoking spatial extent are reviewed. Finally, some effects of early life 
are discussed that contributed to modifying and regulating atmosphere 
and climate of the Earth, and that could have contributed to the highly 
oxidized state of its crust. 



Glossary 

Chiral, achiral and racemic 

A molecule is chiral if its three-dimensional structure is different from its mirror 
image. Such molecules tend to be optically active and turn the polarization 
plane of linearly polarized light in the right or left handed sense. Correspond- 
ingly, they are referred to as D- and L-forms, which stands for dextrorotatory 
and levorotatory molecules. An achiral molecule is mirror symmetric and does 
not have this property. A substance is racemic if it consists of equally many left 
and right handed molecules. A polymer is said to be isotactic if all its elements 
have the same chirality. 



Enantiomers and enantiomeric excess 

Enantiomcrs denote a pair of chiral molecules that have opposite handedness, 
but are otherwise identical. Enantiomeric excess, usually abbreviated as e.e., is 
a normalized measure of the degree by which one handedness dominates over 
the other one. It is defined as the ratio of the difference to the sum of the two 
concentrations, so e.e. is between —1 and +1. 
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Epimerization and racemization 

Epimcrization is the spontaneous change of handedness of one sub-unit in a 
polymer. Racemization indicates the loss of a preferred handedness in a sub- 
stance. 

Catalysis and auto-catalysis 

Catalysts are agents that lower the reaction barrier. A molecule reacts with the 
catalyst, but at the end of the reaction the catalyst emerges unchanged. This is 
called catalysis. In auto-catalysis the catalyst is a target molecule itself, so this 
process leads to exponential amplification of the concentration of this molecule 
by using some substrate. Biological catalysts are referred to as enzymes. 

Nucleotides and nucleic acids 

Nucleotides are the monomers of nucleic acids, e.g. of RNA (ribonucleic acid) or 
DNA (deoxyribonucleic acid). They contain one of four nucleobases (often just 
called bases) that can pair in a specific way. Nucleotides can form polymers and 
their sequence carries genetic information. One speaks about a polycondensa- 
tion reaction instead of polymerization because one water molecule is removed 
in this step. Other nucleotides of interest include peptide nucleic acid or PNA. 
Here the backbone is made of peptides instead of sugar phosphate. 

Peptides and amino acids 

Amino acids are molecules of the general form NH 3 -CHR-COOH, where R 
stands for the rest, which makes the difference between different amino acids. 
For glycine, the simplest amino acid, we have R=H, so two of the bonds on the 
central C atom are the same and the molecule is therefore chiral. For alanine, 
R=CH 3 , so all four bonds on the central C atom are different, so this molecule 
is chiral. A peptide is a polymer generated through a polycondensation reaction 
of amino acids. Peptides are also referred to as proteins. 

Solar constant and albedo 

The solar constant is the total flux of energy from the Sun above the Earth 
atmosphere. Its current value is S = 1.37 kWm~ 2 , but it has been about 30% 
lower when the solar system was young (10 8 yr old, say), so S is not a constant. 
The albedo A is the fraction that is reflected from the surface of the Earth, e.g. 
by clouds and snow and, to a lesser extend, by land masses and oceans. 

Photosynthesis and carbon fixation 

Photosynthesis uses light to reduce CO2 and to produce oxygen either as free 
molecular oxygen or in some other chemical form. This process removes CO2 
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from the atmosphere and produces biomass, which is written in simplistic form 
as (CH20)„. This process is referred to as carbon fixation. 

Life 

A preliminary definition of life involves replication and death, coupled to a 
metabolism that utilizes any sort of available energy. Life is characterized fur- 
ther by natural selection to adapt to environmental changes and to utilize avail- 
able niches. A proper definition of life is difficult given that all life on Earth can 
be traced back to a single common ancestor. Any definition of life may need to 
be adjusted if extraterrestrial or artificial life are discovered. 

Definition of the subject 

Astrobiology is concerned with questions regarding possible origins of life on 
Earth and elsewhere in the Universe. Although there is presently no detection 
of extraterrestrial life, it is generally assumed that life could be wide-spread pro- 
vided certain conditions of habitability are met. A common implicit hypothesis 
in astrobiology is that life can emerge spontaneously once certain environmental 
conditions are met. This implies that there may well have been multiple gene- 
ses, separated only by global extinction events such as major impacts by other 
celestial bodies pQ. 

Four important discoveries can be named that have provided impetus to the 
field of astrobiology. 

1. More than 300 extrasolar planets have been discovered since 1995, provid- 
ing explicit targets for detecting life outside the solar system. 

2. Recent Mars missions have provided evidence for liquid water on the sur- 
face of Mars in past and possibly even present times. This has fostered 
the search for techniques to detect microbial life on Mars. 

3. On Earth the carbon in very old sedimentary rocks dating back 3.8 Gyr 
ago show a consistently lower 13 C to 12 C abundance ratio. This is normally 
indicative of life. This lends support to the notion that life may have been 
present as soon as the Earth surface became hospitable. 

4. The discovery of extremophiles on Earth has considerably extended the 
definition of habitability to include extreme temperatures, pressures and 
pH values, high salinity as well as high radiation levels. This has raised 
hopes of finding life even in our solar system. 

Astrobiology thus comprises several scientific disciplines: astronomy, geol- 
ogy, chemistry, and biology. Therefore, much of the original literature tends to 
appear in journals in the respective fields. We should also mention that there 
are technological attempts in producing artificial life [3J. While this approach 
is not aimed at reproducing the origin of life on Earth, it may still be useful 
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for feeding our imagination in understanding the transition from nonliving to 
living matter. 

1 Introduction 

Since the early days of nonlinear dynamics and non-equilibrium thermodynamics 
it has been clear that one of the ultimate applications of this theory might be 
to facilitate an understanding of the transition from non-living to living matter. 
The main reason is obviously that living systems are very far from equilibrium — 
as indicated by the high degree of order and hence the low entropy of living 
systems relative to their exterior. 

Already in 1952 Turing |3J proposed the idea that chemical reaction-diffusion 
systems might provide a tool for studying biochemical pattern formation which 
have increased our understanding of the laws of nature far from equilibrium, 
where life occurs. This idea was followed up in the late 1960s by Prigogine [H[S] 
who suggested that dissipative structures have great importance in establishing 
a physical description of living matter. A general theory of autocatalytic molec- 
ular evolution was developed in 1971 by Eigen [6] who argued that in a single 
micro-environment only a single handedness can result from a single event. In 
particular the famous chicken and egg problem that occurs in biology at differ- 
ent levels was identified as a Hopf bifurcation. A Hopf bifurcation describes the 
spontaneous emergence of an oscillating solution once some stability threshold 
has been crossed. The mathematics of this is familiar to any physicist, but it re- 
quires that the equations describing the relevant physics are known. In biology 
it is not even clear that the various phenomena can be described by equations. 
A first detailed attempt in this direction was indeed that of Eigen. However, 
the equations governing the emergence of life are only phenomenological ones. 
Nevertheless, these approaches are invaluable in that they help giving the origin 
of life question a mathematical basis. 

One of the earliest anticipated forms of life that is still similar to present life 
is the RNA world [7], whereby simple RNA molecules with functional behavior 
self-reproduces using genetic information encoded either in the same or in other 
participating RNA molecules. Obviously, there are tremendous difficulties given 
that RNA is too complicated a molecule to be synthesized abiologically. A 
significantly simpler molecule is peptide nucleic acid or PNA [5], where the 
backbone consists of peptide instead of sugar phosphate. Nevertheless, the 
difficulty of producing RNA remains. 

There is no firm idea where on Earth such molecular replication may have 
taken place. Frequently discussed scenarios include hydrothermal vent systems 
j9], but also beach scenarios are discussed that are subject to tides leading to 
cyclic changes in concentration |10j as well as to repeated wetting and drying 

An early experiment that contributed significantly to the research into ori- 
gins of life was the Urey-Miller experiment [12] that demonstrated the sponta- 
neous production of amino acids in a reducing atmosphere consisting of H2O 
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vapor, CH4, NH3, and H2 with an energy supply in the form of sparks. More 
recent experiments also allow for the presence of CO2, which now seems un- 
avoidable on the early Earth given that it is continuously replenished through 
outgassing by volcanoes. 

In the following we review some aspects where there has been considerable 
cross-fertilization between astrobiology and nonlinear dynamics. We begin by 
discussing a phenomenon that is believed to have taken place around the time 
of the origin of life, namely the establishment of a definitive handedness of 
biomolecules that is inherent to DNA and RNA (D-form) and to amino acids 
(L-form). Next we discuss constraints on the evolution of hereditary informa- 
tion, and finally review some models that characterize the alteration of the 
terrestrial environment by early life. In addition to any of these physical effects 
there are random fluctuations leading inevitably to local imbalances between 
the concentrations of molecules of D- and L-form. In the following we discuss 
mechanisms that can lead to an exponential amplification of the enantiomeric 
excess. For a recent review of these ideas see Ref. [13] , 

2 Homochirality 

Theories of a chemical origin of life involve polymerization of nucleotides that 
carry and utilize genetic information. Ribonucleotides possess chirality, i.e. they 
are different from their mirror image. All known life forms use ribonucleotides 
of the so-called D-form (right-handed), as opposed to the L-form (left-handed). 
These two molecules are referred to as opposite enantiomers. In most cases 
these different enantiomers are optically active, i.e. they turn the polarization 
plane of linearly polarized light in a right-handed or left-handed sense. 

Any non-enzymatic synthesis of ribonucleotides would have produced a mix- 
ture of equally many right- and left-handed building blocks. Technically this is 
referred to as a racemic mixture of these molecules. However, experimentally 
it is known that in a racemic mixture of mononucleotides the polymerization 
is quickly terminated after the first or second polymerization step [14] . This is 
generally referred to as enantiomeric cross-inhibition, which was long thought 
to be a serious obstacle to a chemical origin of life. It was therefore thought to 
be necessary that life evolved only in a homochiral environment. Moreover, it 
would then be necessary that the degree of enantiomeric purity must have been 
very high. This is important because it rules out a number of physical mech- 
anisms based on the enantioselective effects of circularly polarized radiation, 
magnetic fields, and the parity-breaking property of the electroweak force. 

2.1 The Frank- mechanism 

A general mechanism for producing complete homochirality was proposed in 
1953 by Frank [T5| based on the assumed effects of auto- catalysis and what he 
called mutual antagonism. In fact, the enantiomeric cross-inhibition mentioned 
above can be thought of as a possible example of mutual antagonism. The 
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model of Frank is characterized by the following set of three reactions: 

D + S ^ D + D, (1) 

L + S^L + L, (2) 

D + L^DL, (3) 

where D and L denote monomers of the two enantiomers, S is a substrate 
from which monomers could be formed via auto-catalysis, and DL are inactive 
dimers that are lost from the system. (At this level of simplification no distinc- 
tion between DL and LD is made. This simplification will later be relaxed.) 
The parameters kc and fc/ characterize the reaction speeds. These reactions 
translate to the following set of equations for the concentrations of D, L, S, and 
DL, 

±[D] = +k c [S}[D}-k I [D}[L], (4) 
d [L] = +kc[S}[L)-k I [D}[L], (5) 



dt 
d 
dt 

d 

dt 



[S] = -kc[S]([D] + [L]), (6) 
DL] = +2ki[D][L]. (7) 



These equations imply that the total mass of all building blocks (including the 
substrate), is constant, i.e. [D] + [L] + [S] + [DL] = const = M. 

This system of equations describes the continued autocatalytic production of 
DL, D and L until the substrate S is exhausted, i.e. [S] = 0. However, as long 
as [S] is still finite, the asymmetry, A = [D] — [L], grows quasi-exponentially, 
proportional to exp(/[S'] dt). A numerical example of this is shown in Fig. [TJ 

In the numerical example above we started with very small initial concen- 
trations. Another possibility is to start with a perturbed racemic solution. The 
racemic solution is given by [D] = [L] = X/kj, where A = kc [S] is the instan- 
taneous growth rate due to autocatalysis. Under the assumption that A can 
be treated as a constant (i.e. when the system is still nearly racemic), a linear 
stability analysis shows that the enantiomeric excess, 

ee J^'-W (8) 
^ [D] + [L] W 

grows exponentially. This means that the racemic solution is unstable and that 
the mechanism for achieving homochirality is based on a linear instability. 

2.2 Continued polymerization 

There is a priori no good reason to permit the production of heterochiral dimers 
DL, but not of homochiral dimers DD and LL, i.e. 

D + D DD, (9) 
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Figure 1: Solution of Eqs. (J4j)— ((6j) for hi = kc- Both D and L grow exponen- 
tially until [D] + [L] becomes comparable to the constantly declining substrate 
concentration [S]. At the same time the production of DL removes an equal 
amount of D and L, but this effect affects mostly those enantiomers that are 
already in the minority. In this calculation an initial asymmetry (here 10%) of 
[D] — [L] grows until saturation. At the end, [D] has reached 100% enantiomeric 
excess, but this happened at the expense of producing a large number of inactive 
heterochiral dimers DL. 



L + L^LL. (10) 

The importance of such reactions was stressed in a review by Blackmond [16] . 
who also introduced an additional modification that consists in the assumption 
that not the monomers, but the homochiral dimers DD and LL catalyze the 
production of monomers, i.e. reactions fl} and (|2|) are replaced by 

DD + S ^ DD + D, (11) 
LL + S LL + L. (12) 

This model is similar to the original Frank model provided there is a way of 
getting rid of those homochiral dimers that are in the minority. This requires 
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enantiomeric cross-inhibition for dimers to form heterochiral trimers, i.e. we 
need the additional reactions 

DD + L DDL, (13) 

LL + D LLD. (14) 

A solution to the corresponding reaction equations is given in Fig. [2] Reac- 
tion calorimetry is able to give support to the assumption that dimers and not 
monomers are the relevant catalysts [IB]. This seems to apply in particular to 
the first autocatalytic reaction ever found that enhances enantiomeric excess 
[T?] . In this reaction (sometimes referred to as the Soai reaction) the substrate 
is pyridine-3-carbaldehyde and the chiral molecule of cither D- or L-form is 3- 
pyridyl alkanol which thus acts as an asymmetric autocatalyst to produce more 
of itself. In this reaction, however, dialkylzinc acts as an additional achiral cat- 
alysts. While the Soai reaction is important as a first explicit example of an 
autocatalytic reaction that enhances the enantiomeric excess, it is not normally 
regarded as directly important for astrobiology. 

The polymerization model has been developed by Sandars [18 , who included 
arbitrarily many polymerization steps of the form 



D n +D-=*+D n+ i, (15) 

L n + L-^L n+1 , (16) 

D n + L D n L, (17) 

L n + D^ L n D. (18) 



The basic outcome of this and similar models is always the same as in the 
original Frank model, except that the polymerization model is also capable 
of displaying interesting wave-like dynamics in time-dependent histograms of 
different polymers [15] . 



2.3 Spatially extended models 

In reality there are limits as to how much a system can be considered fully mixed. 
In general, [D] and [L] should be functions of time and space, i.e. [D] = [D](t,x) 
and [L] — [L] (t, x). Assuming that there is only molecular diffusion, the relevant 
reaction equations are to be supplemented by additional diffusion terms, 

^[D n ]=R^+ K V 2 [D n ], (19) 

j t [L n ] =R^+KV 2 [L n ], (20) 

where and R^ L are the right hand sides of the reaction equations. 

If there were only one type of handedness, the resulting equation would be 
reminiscent of the Fisher equation [20] . 

^ = A(1-/)/ + kV 2 /, (21) 
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Figure 2: Solution of Eqs. (j!])-© supplemented by the kinetic equations corre- 
sponding to the reactions Eqs. ©-(O, for k$ = ki = kc- Again, an initial 10% 
asymmetry of [D] — [L] grows until [D] + [L] becomes comparable to the con- 
stantly declining substrate concentration [S]. The monomers polymerize into 
dimers DD and LL. Toward the end, [DD] reaches a 100% enantiomeric excess 

which admits propagating front solutions with front speed Wf ron t — 2 V VcA. Here, 
/ could represent the local concentration of some disease in models of the spread 
of epidemics, for example. 

In the present case there are two fields, one of each handedness. It is instruc- 
tive to refer to these fields as populations which is suggestive of their ability to 
replicate, migrate, become extinct, and to compete against a population of op- 
posite handedness. Each population is able to expand into unpopulated space 
at a speed given approximately by Vf r0 nt, but once two opposing handednesses 
come into contact, there is an impasse and the propagation comes to a halt. A 
snapshot of a one-dimensional model illustrating the polymer length as a func- 
tion of position is shown in Fig. [3] for populations of opposite handedness that 
have come into contact. 

The overall dynamics of symmetry breaking is well characterized by a low 
order truncation, where the model is truncated at n — 2 and the evolution of 
the n = 1 modes is assumed to be enslaved by the evolution of the n = 2 modes 
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Figure 3: Color/gray scale plots of [D n ] and [L n ] for t/Tdie = 0.8 as a function 
of x and n, and the corresponding dependencies of Yl n =i n [-C«] an d 12n=i n [Ln] 
(solid line), compared with ^N^2^ =1 [D n ] and ^^ =1 [L„] (dashed line), and 
±A 2 [L>i] and |iV 2 [Li] (dotted line), all in units of (Q/k s ) 1/2 . The normalized 
diffusivity is k/(L 2 X ) = 10~ 2 and N = 20. Adapted from Ref. [2T] , 



|19j . An example of such a solution is shown in Fig.|4j which shows the evolution 
in a space-time diagram, where two populations of opposite handedness expand 
into unpopulated space until two opposite populations come into contact. 

In two and three dimensions a front between two opposing enantiomers is 
in general curved, in which case it can propagate diffusively in the direction of 
curvature. This is caused by the fact that the inner front between two popula- 
tions is slightly shorter than the outer one. (Only the immediate proximity of a 
front matters; what lies behind it is irrelevant if it is of the same handedness.) 
Indeed, on a two-dimensional surface the inner front is by 2ird shorter than the 
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Figure 4: Profile of [D](x,t*) and space-time diagram of [D](x,t) for the one- 
dimensional problem without advection and an initial perturbation correspond- 
ing to a weak (amplitude 0.01) right-handed excess at x/L = 0.1 (marked in 
white or yellow) and a somewhat stronger (amplitude 0.3) left-handed excess 
at x/L = —0.1 (marked in dark or blue). Note the propagation of fronts with 
constant speed if the exterior is racemic (i.e. [D] = [L] = 1/2, shown in medium 
shades or red) and a nonpropagating front when the chirality is opposite on the 
two sides of the front. The normalized diffusivity is k/(L 2 \q) — 10~ 2 , i.e. the 
same as in Fig. [3] Adapted from Ref. [21] , 



outer one — independent of radius. Here, d is the front thickness, which is of the 
order of d ~ (k/A) 1 / 2 . 

It turns out that in two dimensions the rate of change of the integrated 
asymmetry, A — J([D] — [L]) d 2 x, depends only on the number of topologically 
distinct rings or islands. Once an island is wiped out, the rate of change of A 
changes abruptly and stays then constant until the next island gets wiped out. 
So the enantiomeric excess, 

J([D]-[L])^x 
C -°-- J([D] + [L])d^ (22) 

increases with time in a piecewise linear fashion. 

Even if at each point homochirality could be reached rapidly (time scale 
A -1 ), global homochirality requires that one population wipes out the other one 
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Figure 5: The evolution of molecules of D- and L-forms is shown in the upper 
and lower panels, respectively. Note the tendency toward complete homochi- 
ralization by gradually filling up isolated islands with the chirality of the sur- 
rounding molecules. The dimensionless times are 50, 250, and 1750 from left to 
right. Courtesy of Shibata et al. [22] . 

completely. Diffusion is usually too slow to lead to any significant mixing and 
hence to global homochirality. However, there could be circumstances where 
such mixing is sped up by something like "turbulent" transport. In the case of 
the Earth the slowest relevant transport is in the Earth mantle, part of which is 
now associated with what is called the deep biosphere. Assuming that multiple 
geneses of life is possible, this would raise the question whether a simultaneous 
co-existence of different handednesses on different parts of the early Earth would 
have been possible. It is however unclear whether this possibility could have 
left any traces that would still be detectable today. 

Another approach to solving the problem of spatially extended chemistry 
is by means of cellular automata. In this approach points on a mesh can take 
different states corresponding to molecules of right or left handedness, achiral 
substrate molecules, or even empty states. An example of such a calculation by 
Shibata et al. [32] is shown in Fig. [5] Again, there are patches of populations 
of opposite handedness that grow and wipe each other out such that in the end 
only one handedness survives. 
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2.4 Epimerization 

An interesting alternative to the Frank-type mechanism is a set of reactions 
based primarily on a phenomenon called epimerization, i.e. the spontaneous 
change of handedness in one part of the polymer. This mechanism is important 
in the chemistry of amino acids. Plasson et al. [53] identified four reactions: 
activation, polymerization, epimerization, and depolymerization as necessary 
ingredients that can, under certain conditions, lead to an instability of the 
racemic state with a bifurcation toward full homochirality. They called this 
the APED model, whose reactions can be summarized as follows: 
A: activation: 

L^L*, D^D*, (23) 

P: polymerization: 

L* + L -1+ LL, D* + D-^> DD, (24) 
L* + D^LD, D*+L^DL, (25) 



E: epimerization: 
D: depolymerization: 



LD -U DD, DL LL, (26) 



LL-^L + L, DD D + D. (27) 



This minimal subset of reactions is shown in Fig. [6] 

Compared with the Frank model, a major advantage of the present model 
is that no hypothetical auto-catalysis is required. Indeed, all these reactions 
exist in principle, although it is as yet unclear what kind of manipulations 
on the environment are required to make all these reactions happen. Another 
advantage is that the system is closed, so no inflow or outflow of matter is 
required. The system is maintained away from equilibrium by energy input 
though the activation of amino acids. 

Given that there is neither auto-catalysis nor enantiomeric cross-inhibition, 
one wonders whether the APED model still shares some similarities with Frank's 
original model. Some degree of similarity is immediately seen by writing the 
APED reactions in sequential form in one line, i.e., 

D* + L^> DL^ LL^ L + L, (28) 

L* + D LD -2-> DD D + D. (29) 

This shows that, as long as the reaction rates for epimerization and depolymer- 
ization are not limiting factors, we have essentially the reactions 

D *^ L and L* ^5 D. (30) 

This way of writing these reactions emphasizes the roles of L and D in catalyzing 
the conversion of D* into L and L* into D, respectively. Just like the mechanism 



13 




Figure 6: Representation of the minimal set of reactions necessary for allowing 
the spontaneous transition to homochirality. Adapted from Ref. [24] . 



of mutual antagonism, these reactions disfavor a racemic state, but instead of 
producing unreactive waste, these reactions produce directly one of two possible 
homochiral states. 

In addition, there are reactions of the form 

L* + L LL L + L, (31) 

D* + D -^-> DD D + D. (32) 

These reactions simulate the autocatalytic conversion of L* into L by L and of 
D* into D by D. Again, linear analysis establishes that the racemic state is 
unstable provided a is in the range < a < 1; sec Refs. [23| [24] . 

In conclusion we can say that the homochirality of life-bearing molecules 
might well have originated from the chemical reactions that lead to their for- 
mation. Thus, the hypothetical RNA world may have been born into an en- 
vironment surrounded by homochiral peptides (as described in Section 12. 4p . 
or, alternatively, homochirality may have emerged as a consequence of enan- 
tiomeric cross-inhibition during the first stages of the RNA world (as discussed 
in Section [2~2]). 

In the next section we discuss some issues regarding possible strategies for 
establishing a primitive information-carrying system. This is also based on 
catalysis, but catalysis in the production of other molecules than itself. 
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Figure 7: Number of survivors as a function of the number of generations 
in a numerical example with an error rate of p = 0.01. This initial number 
of survivors is 1000. Note that for a genome length of N = 70 and 71 the 
population dies out after 1800 and 400 generations, respectively. For N = 69 
and less the number of survivors increases exponentially. This is compatible 
with the survival criterion N (ln2)/0.01 = 69.3, derived from Eq. (f33|) . 

3 Establishing hereditary information 

We have so far ignored the fact that polymers can consist of different amino 
acid or nucleotide units, even though they would all have the same handedness. 
Therefore such molecules could in principle carry information. Once such poly- 
mers can replicate, the question arises how to prevent them from getting extinct 
due to errors in the copying process, and instead to compete against parasites. 
It is generally believed that early self-replicating systems had a substantial error 
rate associated with each replication event. A certain small error rate is obvi- 
ously necessary for facilitating Darwinian evolution by natural selection, but it 
must be small enough to prevent extinction. 

Assuming that with each generation a species produces a offspring where the 
length of the genome is N bits, and that the probability for a copying error at 
any position in the genome is p, the necessary condition for long-term survival 
is given by |3 [25] 

pN<\n<r. (33) 

The significance of this formula is illustrated in Fig. [7] with the help of a nu- 
merical example where the selective advantage (i.e. the multiplication factor) 
is chosen to be a — 2, the error rate is p = 0.01, and four different values of 
N between 68 and 71 are used. In this numerical experiment, a new offspring 
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are produced, but with a probability p an error is introduced at each of the ./V 
positions. After this only the intact copies can produce further offspring, and 
so forth. 

According to Eq. Q33p the maximum genome length is, with the parameters 
of our example, (ln2)/0.01 = 69.3. This is compatible with Fig. [7] which shows 
that the dividing line between extinction and long-term survival is between 
N — 69 and 70. For contemporary genomes N is of the order of 10 8 , and p is of 
the order of 10 -8 [25], or below, depending on the efficiency of error-correcting 
mechanisms that are available in contemporary organisms. 

The first replicating systems are likely to have rather high error rates, and no 
correction mechanism, making it virtually impossible to carry sufficient infor- 
mation for building more complex replicators. This difficulty can be removed by 
invoking the concept of hypercycles [27] , whereby the full genetic information 
is carried collectively by several smaller systems (smaller N), each one small 
enough to obey Eq. (j33|) . Mathematically, such a system can be described by 
the following set of reactions [25] : 

h 2 J,, (34) 

h + 2Ii + (35) 

Assuming furthermore that resources arc limited, the total number of molecules, 
M = J^Jlj], is taken to be constant, i.e. U is assumed to be siphoned off from 
the system at a rate <fi that is independent of i. Mathematically, such a system 
can be described by the following set of ordinary differential equations: 

^[J<] = nih] + fei^P-i] ~ 4>{hl (36) 

where 

i I i 

is the factor that keeps the total number of molecules constant. The kinetic 
coefficient ri models the residual effects of birth and death, while ki is the kinetic 
coefficient for the catalytic production of Ii, where acts as a catalyst. The 
evolution of number densities in a model of five hypercycles is shown in Fig. [5J 
for a case where all ki = k and = r are chosen to be the same for all values 
of i. 

An interesting situation arises when the effects of parasites are included. 
Boerlijst & Hogeweg [28] considered an example where a parasite is coupled 
to I2; see Eq. ©. The effect on the above model is shown in Fig. ITD1 where 
fcpara = 2fc and r para = r has been chosen. One sees that not much happens 
for a long time when the parasite is turned on. This is because the parasite 
has to grow to a level where it can affect the entire system. When this point 
is reached, all components of the system decay exponentially — including the 
parasite itself. Unfortunately, the system can never recover from this disaster, 
so the hypercycle theory seems to have a problem. 
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Figure 8: Evolution of the number densities of five hypercycles with equal 
parameters. Note that peaks of I\ (solid line) are followed by peaks of I2 
(dotted line) and -Z3 (dashed line), and so forth. Time is measured in units of 
and concentrations are measured in units of r/k. 




Figure 9: Sketch showing the coupling of several hypercycles together with a 
parasite coupled to species I2. Courtesy of M.C. Boerlijst & P. Hogeweg [25] , 

Again, spatial extent can significantly alter the situation. Using a cellular 
automata approach, Boerlijst & Hogeweg 28] showed that the danger of para- 
sitic catastrophes can be eliminated by allowing the offspring to enter unpolluted 
areas faster than the growth of the parasite allows. Interestingly enough, this 
approach tends to produce spiraling interfaces between different species; see 
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Figure 10: Evolution of the number densities of five hypercycles with equal 
parameters and a parasite where £; para = 2ki and r para = r. Time is measured 
in units of r _1 and concentrations are measured in units of r/k. 




Figure 11: Spatial patterns of the limit cycle showing spiral waves around an 
initial obstacle (left) that was removed at later times (right). Courtesy of M.C. 
Boerlijst & P. Hogeweg [25] . 

Fig. HI 

These equations are in their nature similar to other chemical reaction-diffusion 
equations where several different substances catalyze each others reactions. 
A particularly exciting example is the famous Belousov-Zhabotinsky reaction, 
where malonic acid, CH2(COOH)2 is oxidized in the presence of bromate ions, 
BrO^T. To initiate the reaction, cerium is used as catalyst to donate ions, al- 
though other metal ions are also possible. The color depends on the state of the 
cerium as it changes from Ce 3+ to Ce 4+ or, if iron is used, from Fe 2+ to Fe 3+ . 
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The resulting reactions are of the form [30] 



A + Y -^-> X, (38) 
X + Y P, (39) 
B + X + Z, (40) 

2x g, (4i) 

Z y, (42) 

where X=HBr0 2 , y=Br~, Z=Ce 4+ , A = P=BrC>3~, P and Q are reaction 
products that do not contribute further to the reactions, and fc l7 k^ are known 
rate constants. The reactions above lead to kinetic equations of the form 



d[X] 
dt 



= k, [A] [Y] - k 2 [X] [Y] + k 3 [A] [X] - k 4 [X] 2 , (43) 
d[Y] 



dt 

d[Z] 



k 1 [A][Y]-k 2 [X][Y] + k 5 [Z], (44) 
2k 3 [A][X]-k 6 [Z]. (45) 



dt 

This model of reaction equations is called the Oregonator, which refers to the 
affiliation of the authors at the time of publication [31] . 

If spatial extent is included via diffusion terms, this reaction exhibits in cer- 
tain cases spiral patterns, similar to those in the model of Boerlijst & Hogeweg 
[2"B] . In Fig. rrj] we reproduce the pattern obtained by Zhang et al. [32] for 
a slightly modified model consisting of only two partial differential equations. 
Depending on the value of a certain control parameter in their model spiral 
patterns of different size are being produced. An extensive review of the physics 
of pattern formation in different settings is given by Cross and Hohenberg [33] . 

The connection between pattern formation and the origin of life may seem 
rather remote. However, the equations governing chemical pattern formation 
illustrate some of the critical steps that are expected to play a role in origins 
of life. In particular the fact that different chemical compounds catalyze each 
other in a productive manner is an essential property behind the model proposed 
by Eigen. The Belousov-Zhabotinsky reaction also illustrates the phenomenon 
of autocatalysis where, in the presence of A, the molecule X catalyzes the 
production of more X by using A as a substrate and producing Z as additional 
side product. 

The possibility of self-replication has been demonstrated for simple RNA 
molecules by Spiegelman [31] back in the late 1960s. There are now even ex- 
amples of simple peptide chains that can catalyze the production of each other 
[35| . However, a serious shortcoming of any of the above examples discussed 
here is the fact that there is no possibility of natural selection and hence Dar- 
winian evolution. So, as far as the question of the origin of life is concerned this 
pathway must be considered a dead end. 
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Figure 12: Spiral and ring-like patterns for the modified reaction equations by 
Zhang et al. [35]. On the boundaries a no-flux condition has been adopted, i.e. 
the normal components of all gradients vanish. Courtesy of X.Z. Zhang et al. 
[32]. 

In summary one can say that there are similarities in the mathematics of 
producing homochirality and in establishing hereditary information in the com- 
position of the first replicating polymers. However, in the latter case even less 
is known about the detailed nature of such polymers and their catalytic proper- 
ties. A particularly important aspect is the possibility of spatial extent that can 
substantially modify the behavior of any chemical system. In the present case, 
as shown in Ref . [28] , the possibility of spatial extent is critical for stabilizing the 
system against destruction by parasites. The model also exhibits spiral pattern 
formation that has been at the heart of the early work by Prigogine and others 
in connection with early ideas on biogenesis. 

4 Alteration of the environment by early life 

In this last section we discuss some physics problems within astrobiology that 
illustrate how life, once it has formed, might affect the environment of the early 
Earth and how it led to a planet so markedly different from a planet that does 
not harbor life. 
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4.1 Global energy balance of the Earth 

The young Sun was about 30% fainter than today, and yet the young Earth was 
covered with liquid water and had temperatures higher than nowadays. This 
was caused by the presence of greenhouse gases such as water vapor, carbon 
dioxide, and probably also methane. Life is responsible for reducing CO2 to 
compounds of the form (CH20)„ and similar, and to oxidize various minerals 
or to produce O2. The resulting decrease of CO2 weakens the greenhouse effect, 
so in this sense the emergence of life has essentially a cooling effect on the 
Earth's overall climate. 

Without atmosphere, the planet would cool like a black body at a rate pro- 
portional to the local flux ctsbT 4 , where ctsb is the Stefan-Boltzmann constant. 
Integrated over the entire surface of the planet, this corresponds to a loss of 
47ri?|cr S BT 4 , which would need to be balanced against the rate of energy re- 
ceived by solar radiation. The solar "constant" is S = 1.37kWm~ 2 and the 
total energy projected onto the disk of the Earth is (1 — A)ttR^S, where A is 
the albedo, i.e. the fraction of energy reflected from the Earth. The resulting 
blackbody temperature would be 

1/4 

(46) 

Using A = 0.3 and a SB = 5.67 x lO^Wm" 2 ^ 1 ), the temperature of the 
Earth would be 255 K or about —18 C. 

In the presence of an atmosphere the rate of cooling is modified to <7sbT* s , 
where T c ff is the effective temperature equivalent to that of a black body. A 
positive greenhouse effect corresponds to T c s < T, so the cooling is reduced and 
the atmosphere heats up according to the vertically integrated energy equation 

C^ = (l-A)|-a SB T c 4 ff , (47) 

where C is the vertically integrated specific heat. 

The value of the effective temperature can be obtained from a radiative 



Tn = 



(1-A) 



S 



4cr SE 
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transfer calculation. A simplified model calculation^ yields 




(48) 



where To is the surface temperature, £ is an averaged mean free path of photons 
and ^crit is the critical value above which there is a positive greenhouse effect. 
Again, a simplified calculation suggests ^ C rit = 3JJ/16 w 0.197J, where H = 
1ZT / (fig) ~ 8 km is the pressure scale height of the atmosphere. So, an increase 
in opacity leads to a decrease of the cooling and hence to an increase in the 
surface temperature. 

Another interpretation is to say that the greenhouse gases shift the radiat- 
ing surface by a certain amount, £ s , upward. The value of £ s is related to I, 
Ditlevsen [36] uses £ g — 3 km. He also noted that a more accurate lapse rate of 
the temperature is dT/dz = lOKkm -1 instead of Tq/H = 40Kkm~ 1 , so that 
the temperature gain caused by greenhouse gases is £ g x dT/dz = 30 K. The 
reason for a shallower temperature gradient is the presence of convection that 
causes the specific entropy s to be nearly constant with height. In that case 
the temperature gradient is just the adiabatic one, (dT/dz) a d = g/c p , where c p 
is the specific heat at fixed pressure, which in turn is related to the universal 
gas constant and the specific weight via TZ/ /i = c p — c v , where c v is the specific 
heat at fixed volume and c p [c v = 7 is the ratio of specific heats. With these 
formulae one does indeed gecl 



where we have used 7 = 7/5 for air molecules with 5 degrees of freedom (3 for 
translation and 2 for rotation). 

At certain times over the history of the Earth other greenhouse gases such as 
methane may have played an important role in keeping the Earth above freezing 

1 Under the assumption of local isotropy (Eddington approximation) radiative equilibrium 
implies that the flux is proportional to the negative gradient of the radiative energy density 
aT 4 , where a is the radiation-density constant, so 



Here, c is the speed of light and I is the mean free path of photons. The latter can be expressed 
in terms of the opacity k and the density p as I = (ftp) - . Hydrostatic equilibrium can be 
written in the form 



where g is the gravitational acceleration, 1Z is the universal gas constant, and p is the mean 
molecular weight. These equations can be solved by a polytrope, i.e. T = To(l — z/H) and 
P = Po(l — Z /H) 3 , where 2 is the distance from the surface and H is the vertical pressure 
scale height. This leads to a condition of the form Eq. I|48|l where i cr it = 3H/16 ~ 0.19-ff is 
the critical mean free path of photons. 

2 Hydrostatic equilibrium can be written as — p~ x Vp — V0 = 0, where p is the pres- 
sure and <j> = gz + const is the gravitational potential. Using the thermodynamic relation 
— Vp = -V/i+TVs = 0, where h = c p T is the specific enthalpy and Vs = for adiabatic 
stratification, we have d(c p T)/dz = g. 




(49) 



F = — 



i C £V(aT 4 ). 
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temperatures. Indeed, the burial of oxides in the crust allowed methane to build 
up in the atmosphere, which may have led to concentrations of a few thousand 
times greater than modern levels. UV radiation in the upper atmosphere breaks 
up methane into its components, letting H2 to escape into space, leading to a 
net gain of oxygen, that comes ultimately from H2O. 

According to a model of Catling et al. [37] methane (CH4) may have been 
important 2. 7. ..2. 3 Gyr ago, just before the famous Snowball Earth deep freeze of 
the Earth [35]. As discussed above, the associated loss of hydrogen may have led 
to a gradual accumulation of oxygen in the atmosphere, which then terminated 
the methane era and led to the Snowball Earth event. This event lasted until 
the continuous CO2 production from volcanoes accumulated to large amounts so 
that the resulting greenhouse effect became sufficient to initiate partial melting 
of the ice cover. 

4.2 Response to changes in greenhouse gases 

As was known from global climate models [13] and later from simplified models 
[H] using Eq. ([4"T]) with a relatively simple piecewise linear temperature depen- 
dence of A(T), there can be three different equilibrium temperatures. This is 
illustrated in Fig. [T3] where we compare the graph of crsBT c 4 ff versus surface 
temperature To with the net radiation (1 — A)S/4. Here, A(T) has been ar- 
ranged such that A = Ah ot for T > T max (corresponding to no ice coverage) 
A = A cold for T < T min (corresponding to full ice coverage). 

Ditlevsen [35] used Eq. 1[17J) to study the response of the system to variable 
greenhouse gas concentrations. As the amount of CO2 increases, the equilibrium 
temperature increases. Obviously, when the system was on the lower fixed point 
initially, there must be a critical CO2 concentration above which the solution 
will jump discontinuously to the upper branch; see Fig. 1141 

It is generally accepted that the rate of weathering increases with increasing 
temperature. This provides a stabilizing effect on the climate. As T increases, 
the rate of weathering increases, removing more CO2 from the atmosphere, 
reducing the greenhouse effect, and thus leading to cooling. Ditlevsen [36] in- 
troduced the assumption that there is a continuous source of CO2 through out- 
gassing from volcanoes and a temperature-dependent sink of CO2 from weath- 
ering when T exceeds a critical temperature T w , but no weathering for T < T w 
[4"S"] . This leads to a self-regulating effect for T < T w , which Ditlevsen calls a 
greenhouse thermostat. Whenever T < T w , since there is then no weathering 
and hence no sink of CO2, greenhouse gases will build up until the Earth's tem- 
perature has reached the value T w ; see Fig. [15] This is the mechanism that is 
believed to have caused the early Earth to be above freezing through most of its 
history — with the exception of intermediate Snowball Earth-like events that are 
caused by the emergence of other sinks of greenhouse gases such as the onset of 
aerobic photosynthesis or the enhanced formation of mountain topography that 
leads to an increase in the erosion rate and hence the weathering. 
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Figure 13: Plot of asBT^g versus surface temperature To for three different 
greenhouse temperature shifts, T g , compared with the net radiation i(l — A) 
for a simple piecewise linear function A(T). Courtesy of P. Ditlevsen [36] , 

4.3 The Daisyworld model 

Life can also affect the planet's albedo, as has been demonstrated by Love- 
lock [39] in his Daisyworld model. For a tutorial on the Daisyworld model see 
Ref. [40]. This model also makes use of Eq. (j47|) . but now the planet's albedo A 
is affected by the plant population which is simplistically represented by black 
and white plants or flowers (daisies) with local albedos A\ and A2, respectively. 
So the total albedo is a weighted average of the form 

3 

A = ^a l A l , (50) 

i=l 

where A3 is the albedo of the unpopulated surface. The weights 04 depend 
on the surface coverage of the respective regions and obey evolution equations 
that are in turn governed by a temperature-dependent growth term, /3(T), and 
a fixed death rate, 7, so the resulting equations for the rate of change of the 
albedo are 

n cv ■ 

^ = [/3(T)- 7 H, (51) 

where i = 1 for black and i — 2 for white plant populations, (3(T) is assumed 
to be different from zero in the range T m \ n < T < T max with a maximum at 
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Figure 14: Equilibrium temperature as a function of CO2 concentration. Cour- 
tesy of P. Ditlevsen [36] . 

?aver = 5 (T" m i n + T max ). The weight for the unpopulated surface follows from 
the normalization Y] a,j = 1, so 0:3 = 1 — a.\ — a-i. 

The temperatures are higher in the regions of black plants and lower in 
regions of white plants according to the formula 

T^(A-A t )q + T \ (52) 

where q is a parameter that must be smaller than a critical value, 

q < gcrit = S/(4ohb), (53) 

in order that heat flows against the temperature gradient ,42 . For q = q cl it 
the temperature is uniform for different values of A4, while for q < q cr it the 
regions of high albedo are cooler and those of low albedo warmer. Note also 
that Eq. (fB"2")) preserves heat balance, i.e. Y2 a i^i = ^o- 

The important point in the Daisy world model is the fact that, for a certain 
range of S, the surface temperature of the planet, T , is stabilized to be in a 
certain temperature range around the optimal value close to T avGr ; see Fig. [T5] 

Saunders makes another remarkable point. He showed that by changing the 
model to allow for Darwinian evolution in such a way that each plant species 
works with an optimized temperature dependence, so /3(Tj) — > Pi(Ti) is modified 
to become dependent on i, the overall result changes only very little. More im- 
portantly, the range over which the model can stabilize the planet's temperature 
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Figure 15: Dependence of surface temperature on time under the assumption 
of a continuous source of outgassing of CO2 and the onset of a CO2 sink for 
T>T W . Courtesy of P. Ditlevsen [36]. 

shrinks, making the planet as a whole more vulnerable. Although the amount of 
shrinkage is small, it emphasizes the dangers associated with adopting changes 
that only lead to short-term benefits. Saunders emphasizes in his work that the 
ability of life to regulate the surface temperature of a planet is not associated 
with natural selection as in the concept of Darwinian evolution. More generally 
he warns therefore that not everything that is to an advantage needs to be the 
result of natural selection [42] . 

4.4 Oxidation of the Earth's crust 

It has recently been proposed that, in addition to the effects discussed above, 
life may have profound effects also on the Earth's crust. A possible scenario of 
such a suggestion has recently been discussed by Rosing et al. [IB] . The idea is 
that photosynthetic life may tap large amounts of solar energy that were used 
to reduce carbon from CO2 to compounds of the form (CI^O)^ and similar, 
via reactions of the form 

C0 2 +H 2 + /w;^ CH2O + O2, (54) 

where hv denotes the energy taken from solar radiation. Furthermore, and 
even more surprisingly, the oxygen produced by photosynthesis may have been 
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Figure 16: Temperature (in Celsius) versus relative irradiation (normalized to 
the average temperature). Note that the temperature is stabilized around the 
value T aver , provided the energy input S is within a certain range. Courtesy of 
P.T. Saunders @2|. 

critical in oxidizing iron in the continental crust. Although other factors also 
plaid a substantial role, it is clear that biological processes can speed up the 
oxidation process substantially. Comparing the oceanic crust with continental 
crust, a major difference is the enhanced fraction of Si02 (57% in the continental 
crust compared to 50% in the oceanic crust). 

With granite being one of the lightest rock types, it was eventually able to 
escape subduction and to produce stable continents about 3.8 Gyr ago. This 
is also the time of the oldest rock findings on Earth. Given that the rise of 
continents on the early Earth is associated with granite formation, the presence 
of granite on silicon-bearing rocky planets might thus be a possible biomarker 
for photosynthesis @B]. 

Although this idea is speculative, it may be supported quantitatively as 
follows. Firstly, the present day production rate of organically fixated carbon is 
estimated to be 9 x 10 15 molCyr _1 [4"Tll4"8] . The amount of energy required for 
this can be calculated by using the fact that it costs 477 kJ to transfer one mol 
carbon to hexose. The energy required for this is then 300 mWm -2 . Rosing et 
al. [46] argue that this amount could be supplied by only 0.1% of the effective 
solar energy flux, S/4. Assuming that the amount of carbon burial, relevant to 
estimating the usable fraction of oxygen for iron oxidation, is also about 0.1%, 
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this corresponds to about 10 13 mol Cyr -1 . This would yield a comparable iron 
oxidation rate. Rosing et al. 46J argue further that the annual basalt production 
contributes about 10 14 mol Feyr -1 , so a fraction of the magmatic iron flux could 
be used for building up the mantle reservoir of ferric iron. 

In conclusion, the presence of life can lead to significant alterations of the 
planet in a number of different ways, as is quite clearly demonstrated by some of 
the differences between Earth and its neighboring planets Venus and Mars. Only 
the Earth has extensive reservoirs of oxygen and of granite. Within limits, the 
presence of life on a planet can also have a stabilizing effect on its climate. The 
relevant mathematical modeling of some of these processes resembles in many 
ways those encountered earlier in studies of homochirality and of the spread of 
hereditary information on the early Earth. 

5 Conclusions 

Astrobiology has developed into a rapidly growing research field involving exper- 
tise from a number of neighboring disciplines. Nonlinear dynamics and noncqui- 
librium thermodynamics find applications in all these subfields. Here, we have 
elaborated on a few such aspects. Closest to the onset of life is perhaps the 
emergence of homochirality of biomolecules. Given that RNA has been proven 
to form longer polymers only in a homochiral environment, one would expect 
that homochirality must be a prerequisite to the emergence of life at the level of 
a replicating RNA world. On the other hand, the very mechanism causing the 
polymerization to terminate, namely the enantiomeric cross-inhibition, can also 
be the mechanism responsible for causing 100% homochirality by destroying 
RNA molecules whose chirality is already in the minority. This would however 
requite the possibility of auto-catalysis, which can be avoided in another sce- 
nario where a closed peptide system is kept away from equilibrium by continuous 
activation of amino acids. 

Chemically speaking, the stabilization of a definite chirality is in some models 
similar to the subsequent establishment of hereditary information in that catal- 
ysis plays a crucial role. Furthermore, in both cases the possibility of chemistry 
in an extended system is crucial. On the one hand, spatial extent gives rise to 
the possibility of coexistence of life forms of opposite handedness on the early 
Earth. On the other hand, spatial extent can be critical in allowing the system 
to find unpopulated locations faster than being overwhelmed by the effects of 
parasites that tap the same resources that are required for the maintenance and 
development of hereditary information. 

Finally, life is invariably coupled to come kind of metabolism that is ulti- 
mately powered by solar energy. This clearly affects the environment by reduc- 
ing carbon and oxidizing the crust of the Earth and, over the last two billion 
years also the atmosphere. How much these alterations of the environment are 
due to biological processes is less obvious. However, it is clear that biological 
factors greatly speed up weathering on the Earth. The extent of biologically 
induced alterations of the continental crust, for example, may therefore best be 
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tested using quantitatively accurate model calculations. The outcome may ulti- 
mately hinge on energetic considerations and on the efficiency of photosynthesis 
as a solar energy collector. 

With the scope of being able to explore in the near future not only the 
planets and other celestial bodies in the solar system in much more detail, but 
also planets of other planetary systems, the research in astrobiology quickly 
develops into a field that will be driven more and more by new data, making 
this field less susceptible to speculation. It is therefore important to be prepared 
for upcoming discoveries in this field. Finally, it should be emphasized that 
astrobiology is efficient in communicating science to the general public, which 
may provide additional boost to the field. 
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